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Abstract 



We consider a S* = 1 kagome Ising model with triquadratic interactions 
around each triangular face of the kagome lattice, single-ion anisotropy and 
an applied magnetic field. A mapping establishes an equivalence between the 
magnetic canonical partition function of the model and the grand canonical 
partition function of a kagome lattice-gas model with localized three-particle 
interactions. Since exact phase diagrams are known for condensation in the 
one-parameter lattice-gas model, the mapping directly provides the corre- 
sponding exact phase diagrams of the three-parameter S = 1 Ising model. As 
anisotropy competes with interactions, results include the appearance of con- 
fluent singularities effecting changes in the topology of the phase diagrams, 
phase boundary curves (magnetic field vs temperature) with purely positive 
or negative slopes as well as intermediate cases showing nonmonotonicity, and 
coexistence curves (magnetization vs temperature) with varying shapes and 
orientations, in some instances entrapping a homogeneous phase. 



PACS numbers: 05.50+q, 75.10.-b 



I. Introduction 

The S = 1 Ising model continues to be actively investigated with applications in phase 
transitions, critical and multicritical phenomena. The special behaviors of the model are 
largely attributable to the presence of single-ion anisotropy and biquadratic interaction 
terms in the Hamiltonian. A familiar version is the Blume-Emery-Griffiths (BEG) model 
lU having bilinear and biquadratic nearest-neighbor pair interactions and single-ion-type 
uniaxial anisotropy. The three parameter BEG model reduces to an earlier Blume-Capel 
model if the biquadratic interactions are neglected. The BEG model has attracted 
considerable attention since it was originally proposed to describe phase separation and 
superfluid ordering in He^-He^ mixtures exhibiting tricritical behavior when the anisotropy 
is sufficiently strong and competes with the interactions. Subsequently, the model has been 
employed to explore a variety of other phase transition problems 0. Diverse theoretical 
approaches have been adopted 0, but invariably due to severe mathematical difficulties, 
approximation schemes are used in the calculations. 

The free energy (characteristic function) of a lattice-statistical model can be found in 
principle by evaluating the partition function ("sum over states"), a formidable mathemat- 
ical task in a highly cooperative system with macroscopically many degrees of freedom. 
The phase diagrams of the system indicate the locations and other characteristics of the 
mathematical singularities in the thermodynamic (large lattice) limit of the free energy per 
site, /. Approaching the singularities, one examines whether the ordering parameter (ap- 
propriate first-order derivatives of /) vanishes continuously or jumps discontinuously as a 
function of the temperature. The phase transition is then called continuous or discontinuous, 
respectively, and is so signified upon the phase diagrams. 

In general, the phase diagrams of the S = 1 BEG model are not known exactly. However, 
neglecting its bilinear interactions, Griffiths |^ has pointed out that the partition function 
of the remanent two-parameter 5* = 1 Ising model is equivalent to the partition function of 
a standard S = 1/2 Ising model in a field, whose exact phase diagrams are well-known in 
the two-dimensional (d=2) ferromagnetic case. Wu extended the investigations by con- 
sidering a three-parameter 5" = 1 Ising model having a magnetic field, single-ion anisotropy 
and biquadratic interactions, and proceeded to obtain exact phase diagrams of the model 
for ferromagnetic interactions on planar lattices. In the present paper, a S = 1 kagome Ising 
model is taken to have a magnetic field, single-ion anisotropy and triquadratic interactions 
around each triangular face of the kagome lattice. A mapping establishes an equivalence 
between the magnetic canonical partition function of the model and the grand canonical 
partition function of a kagome lattice-gas model with localized triplet interactions. Since 
exact phase diagrams are known for condensation in the one-parameter lattice gas model, the 
mapping directly provides the corresponding exact phase diagrams of the three-parameter 
S = 1 Ising model. As anisotropy competes with interactions, special features of the find- 
ings include the appearance of confiuent singularities effecting changes in the topology of the 
phase diagrams, phase boundary curves (magnetic field vs temperature) with purely positive 
or negative slopes as well as intermediate cases showing nonmonotonicity, and coexistence 
curves (magnetization vs temperature) with varying shapes and orientations, in some in- 
stances entrapping a homogeneous phase. To our knowledge, the exact phase diagrams are 
the first found for any S = 1 planar Ising model with multi-spin interactions. 

As added incentive for these studies, it may be instructive to briefiy (albeit heuristically) 
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discuss multi-particle cyclic exchange processes. Letting S/ be a spin operator localized on 
a lattice site /, an isotropic higher-order Heisenberg exchange operator 

(S,-S,)(S,-Sfc)(S,-S,) (1.1) 

can be associated with cyclic exchange around the vertex sites i, j, of a triangle. Although 
not having the precise electron spin operator forms (1.1), analogous multi-particle exchange 
models may be relevant for particles with hard core interactions where direct exchange 
is hindered at high densities with the inference that the simplest exchange involves three 
particles in a ring following one another in a circular motion. An analogy can be experienced 
in an overcrowded bus or elevator where it is easier for a few neighbors to revolve in unison 
than it is for two neighbors to exchange their positions. Three-particle cyclic exchanges 
are known to be important, e.g., in quantum solids and liquids like helium 0. The highly 
anisotropic or Ising version of (1.1) becomes 

{s!s]){s]s^,msn = {s:ns]Y{s^,f, (1.2) 

which is the operator form of the S = 1 triquadratic Ising interactions currently consid- 
ered. The longitudinal Ising form (1.2) is essentially the purely potential energy part of 
the Heisenberg form (1.1). Although Ising model descriptions may not be fully realistic, 
they can contain important vestiges of the underlying complex problem, and offer valuable 
insights amid these complexities. Exact solutions in the simplified models can illustrate the 
role of multi-spin interactions in general and can in particular highlight any possible novel 
effects of the multi-spin interactions that may not be accessible or clearly delineated using, 
e.g., perturbation theory, closed- form approximations, or finite-size numerical simulations. 

The paper is organized as follows. Section II presents the S = 1 kagome Ising model 
whose partition function is suitably transformed using idempotent lattice-gas variables. In 
Section III, the central mapping concepts in the theory are established. Applying the map- 
ping relations, phase boundary curves of the present S = 1 Ising model are determined 
in Section IV as are the companion coexistence curves in Section V. Section VI contains 
concluding comments on the theory and results. 

II. Model Hamiltonian and partition function 

The kagome lattice (Japanese woven bamboo pattern) is a two-dimensional periodic 
array of equilateral triangles and regular hexagons (Fig. 1) thus also called the 3-6 lattice. 
The total numbers of triangles and hexagons are in a 2:1 ratio, with the corner- sharing 
triangles having dual orientations, say, up- or down-ward pointing. The lattice is regular 
(all sites equivalent, all bonds equivalent) and may be termed "close packed" since it contains 
elementary polygons having an odd number of sides, viz., triangles. Note that the kagome 
lattice has the same coordination number 4 as the square lattice, the latter being "loose 
packed" . 

Consider the following 5* = 1 Ising (dimensionless) Hamiltonian defined on the planar 
kagome lattice of Af sites 

-(3n = hJ2s, + DY.s^ + Q, Yl s^s]sl (2.1) 

i i <i,j,k> 

where Si = 0,±1, I = 1, ■■■ ,Af are site-localized spin-1 Ising variables, the summation 
is taken over all lattice sites, J2<i,j,k> is over all triplets of sites belonging to elementary 
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triangles, and (3 = l/ksT with being the Boltzmann constant and T the absolute tem- 
perature. Also, h = (3jlhz is a uniform (dimensionless) magnetic field with jl being the 
electronic intrinsic magnetic moment, D = (31) a single- ion-type uniaxial (dimensionless) 
anisotropy parameter, and Qs = /^Js > a three-spin (dimensionless) interaction parame- 
ter. In (2.1), the "symmetry breaking" impressed field is longitudinal (along z-axis) with 
the transverse x — y plane taken as the plane of the two-dimensional lattice. The magnetic 
field hz and the anisotropy field P couple to the dipole and quadrupole moments of the 
spin assembly, respectively. In the paper, the lattice-statistical model (2.1) is termed a S=l 
triquadratic Ising model. The primary purpose of the theoretical investigations is to deduce 
exact phase diagrams of the model by employing transformation and mapping techniques 
upon its partition function. 

The magnetic canonical partition function >Z(/i, Qs) of the spin system (2.1) is defined 

as 

Z{h,D,Q;) = Y.e-^'^ (2.2) 

where the summation is taken over all possible values of the set {Si} of spin variables. 
Explicitly entering (2.1), the partition function (2.2) may be written as 

= ^(2cosh/i)I^«"' X e^^«"'+'53E<,,,,.>"«"."fc 

{rii} 

{Ui} 

having made repeated use of the partial trace identity 

E e'''fiSf)= E (2cosh/.)"7(n.), (2.4) 

5^=0,^11 ni=0,l 

which is readily established by identifying the terms corresponding to Si = ±1(0) on the 
LHS with those corresponding to rii = 1(0) on the RHS. The enlistment of idempotent 
variables rij, z = 1, ■ ■ ■ , A/" in (2.3) will be useful to directly establish a mapping between the 
partition functions of the S = 1 triquadratic Ising model and a triplet-interaction kagome 
lattice gas model whose phase diagrams for condensation are known exactly. 
III. Lattice gas representation of S" = 1 triquadratic Ising model 
Consider a lattice gas of N atoms upon the kagome lattice of J\f sites with the (dimen- 
sionless) Hamiltonian 

- PHig = Ks J2 ^i^j^k, (3.1) 

where ii's = /3e3 with £3 > being the strength parameter of the short-range attractive 
triplet interaction, the sum is taken over all elementary triangles, and the idempotent site- 
occupation numbers are defined as 

_ J 1, site / occupied , s 

\ 0, site / empty. 
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In (3.1), an infinitely-strong (hard core) repulsive pair potential has also been tacitly assumed 
for atoms on the same site, thereby preventing multiple occupancy of any site as reflected 
in the occupation numbers. 

In the usual context of the grand canonical ensemble, we introduce 

}1 = 1-Lig- liN (3.3) 

where is the chemical potential with N being the conjugate total number of particles 

N = Y.n, . (3.4) 

i 

Using (3.1), (3.3) and (3.4), the grand canonical partition function E{fi,Af,T) is given by 
E{fx,X,T) = J2 e-^"^ = E e^'^^«"'+'^^^<».^.fc>"'"^"^ (3.5) 

{rii} {rii} 

Comparing (2.3) and (3.5), we conclude that 

Z{h,D,Q^) = E{fi,Af,T) (3.6a) 

provided that 

ln(2 cosh h)+D = Q3 = K3. (3.6b) 

The phase diagrams of the triplet-interaction kagome lattice gas (3.1) are known exactly 
0,more definitely, its liquid- vapor phase boundary (chemical potential or pressure versus 
temperature), coexistence curve (density versus temperature) and various critical properties. 
Hence, as shown shortly, the mapping (3.6) affords a convenient method for directly deter- 
mining the corresponding magnetic phase diagrams of the S = 1 triquadratic Ising model 
(2.1). 

IV. Phase boundaries of the S = 1 triquadratic Ising model 

Considering the kagome lattice gas (3.1) with purely three-particle interactions, its liquid- 
vapor phase boundary curve is given by 

/i/es = -Kg"' ln[(e^3 - 1)^/=^ - 1], < K,JK, < 1, (4.1) 

with fi/es being a reduced chemical potential and K^c/K^{= T/Tc) a reduced temperature 
where K^c{'= ea/^s^c) = ln[(2 + v^)^ + 1] = 3.96992 .... The curvilinear phase boundary 
begins at zero temperature with yu/eg = —2/3 and ends (analytically) at a critical point 
whose coordinates are K^c/K^ = yu/cs = /^c/cs = —0.64469 . ... At zero temperature, 
the phase boundary curve fi/es vs T/T^ has a zero slope in accordance with the Clausius- 
Clapeyron equation and third law of thermodynamics. Otherwise, its slope is positive, which 
is more discernible at temperatures closely below the critical temperature. 

To similarly investigate the 5* = 1 triquadratic Ising model (2.1), one directly applies the 
correspondence relation (3.6b) to expression (4.1) thereby yielding 



.gQ3c/x _ 1)2/3 _ 1 2 cosh [^]\, < X < 1, (4.2a) 
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where 

^ = h/Qs{=flK/Jsy, a = D/Q^{^V/J^)- x = Q^JQ^{^T /%)■ 

Qsc = Js/keTc = ln[(2 + Vsf + 1] = 3.96992 • • • . (4.2b) 

Regarding notations, the initial parameters h,D,Q-i in (2.1) are dimensionless through use 
of a thermal energy fc^T scale. Here, for graphical preferences, the associated parameters 
C,C(,x in (4.2) are reduced using an interaction energy J3 scale. 

Expression (4.2) is plotted in Figure 2 as a phase boundary surface x — x(^, a), beginning 
at zero temperature (,x = 0) as a wedge-shaped locus and terminating along a critical {x = 1) 
curve. The phase boundary surface has a number of unusual characteristics, leading to 
several novel features in the phase boundaries and coexistence curves. In the following, we 
will discuss the projections of the surface onto the ^ — ck and $,—x planes which will illustrate 
some of these special features. 

In Figure 3, isotherms of (4.2) in the range < a; < 1 are plotted in ^ — a space ( "field 
space"). The resulting transition region may be viewed geometrically as the projection of 
the phase boundary surface (Figure 2) 'm.^ — a — x space upon the ^ — a plane. In particular, 
(4.2) shows, at ^ = 0, that the zero-temperature {x — 0) isotherm in Figure 3 has its angular 
apex at a^"^ = — 2/3 and that the critical {x — 1) isotherm exhibits a symmetric rounded 
maximum at a;!'"^ = -Qg^^ ln{2[(e'53- - 1)2/^ - 1]} = -0.81929 ■ ■ ■ . Observing Figure 3, it is 
evident that the model only admits phase transitions for anisotropy sufficiently competing 
with interactions {T> < — (2/3)J3). The existence of surface undulations in Figure 2 leads 
projectively to the existence of crossing points among isotherms in Figure 3. For instance, 
the zero temperature {x — 0) and critical temperature {x — 1) isotherms intersect one 
another at {Across, Oicross) — (±0.30171 • • • , —0.96837 • • •) as shown in Appendix A. Note that 
not all isotherms cross at a single point, as shown in Figure 4 by magnifying the boxed region 
of Figure 3. It is also clear from Figure 3 that for cases a^"^ < a < a^"^, only non-critical 
transitions exist, while for cases a > a'^'^^, the model is devoid of any phase transition. 

In practice, a magnetic field parameter ^ is experimentally adjustable, contrasting the 
intrinsic crystal-field anisotropy parameter a. This makes the projection of Figure 2 onto 
the ^ — x plane more relevant from an experimental point of view. As an example, for a given 
specimen, say a = —1.5, and for ^ > 0, the relation (4.2) shows that the phase transitions 
only occur within the range 0.8333 • • • < ^ < 0.8550 • • •, where the lower and upper end 
points of the interval correspond to the zero and critical temperatures, respectively. Since 
the interval does not contain vanishing values of ^, the associated phase transitions for the 
specimen are induced by the applied magnetic field hz- For this example, the relation (4.2) 
is plotted in Figure 5, illustrating the phase boundary (solid curve) to be curvilinear with 
a positive slope at non-zero temperatures. The phase boundary separates two homogeneous 
phases, viz., the longitudinal upper phase associated with the Ising spins preferably aligning 
parallel to the applied field (along the positive 2;-axis) and the transverse lower phase 
associated with the spins preferably lying in the x — y plane of the lattice. The phase 
boundary is the locus of discontinuous transitions from one phase to the other, and the 
terminating point (solid circle) is a critical point associated with a continuous transition. 
In Figure 5, the phase boundary is analytic at the critical point. In fact, the present phase 
boundary curves (magnetic field vs temperature) are analytic at their critical points for 
all a < a™"^. As explained shortly, the phase boundary curve for a — a^"'^ possesses a 
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confluent singularity at criticality (a; = 1). 

The monotonically increasing hehavior of the phase boundary curve in Figure 5, however, 
is not reahzed for all values of the single- ion anisotropy parameter a. Figure 6 shows 
the monotonically increasing phase boundary curve for o; ^ —1 gradually crossing over to 
a monotonically decreasing phase boundary curve for a > o; ~ —.93724. In particular, 
selecting the value across — —0.96837 • • • from within this crossover range, the resulting 
phase boundary curve in Figure 7 exhibits non-monotonic behavior, viz., a shallow rounded 
minimum having coordinates {xmim Cmm) = (0.80577 • • • , 0.30020 • • •). This non-monotonic 
behavior of the phase boundary curves is associated with the region of crossing points in 
Figure 3 (see Appendix A). 

Upon further increasing the values of a, the corresponding ^-coordinates of the critical 
points (solid circles) descend in value, as shown in Figure 8, until the parameter value 
a — aJT"^ = —0.81929 • • •, when the phase boundary curve vanishes as A^e^/"^ (algebraic 
branch point singularity). Here, e = {T^ — T)/Tc is the fractional deviation of temperature 
from its critical value and the calculated amphtude ^4^ = 0.20963 • • • . This singular behavior 
arises as the critical point (^ > 0) shown in Figure 8 meets its reflective image critical 
point < 0) at a confluence point (x, ^) = (1, 0), thereby changing the topological nature 
of the phase diagram. Continuing to increase a beyond its aJJ^"^ value has the effect of 
further shortening the length of the phase boundary curve (comprising solely non-critical 
transitions) whose confluent singularity moves toward the origin, carrying the same exponent 
1/2 but diminishing amplitude. Finally, as expected, the phase boundary curve no longer 
exists for a > —2/3. 

V. Coexistence curves of the S* = 1 triquadratic Ising model 

In the previous section, mapping relations were directly employed to facilitate the finding 
of exact phase boundary curves of the S — 1 triquadratic Ising model. In the present section, 
the companion coexistence curves of the model are determined via composite differentiation 
of the same mapping relations (3.6). More particularly, the (dimensionless) magnetization 
m is represented as 

oh 

_ , 91nS(/x,Ar,T) d{f3fi) 

d{P^) ' dh 
a[ln(2cosh/i) + D] 

" ^ ^ dh 

= p ■ tanh h (5.1) 

where p =< Ui > is the average particle number density oi the triplet interaction kagome lat- 
tice gas (3.1). In the context of coexistence phase diagrams, the correspondence relationship 
(5.1) specializes to 



Kr = Prr(tanh/i)^''^ (5.2) 



where m^'^^ is the longitudinal-transverse coexistence surface in m — a — x space of the 
S — 1 triquadratic Ising model, pf^^ is the liquid-vapor coexistence curve in p — x space 
of the triplet interaction kagome lattice gas and (tanh/i)^''* is evaluated upon the earlier 
determined phase boundary surface (pbs) (4.2), the latter yielding 
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1^1 

1)2/3-1 I , 0<a;<l. (5.3) 



In the product representation (5.2), the a-dependence of m™^^ resides solely in the factor 
{tanh hy^^ as explicitly given in (5.3). For a fixed value of the anisotropy parameter a, the 
expression (5.3) evaluates (tanh/i)^^* along the associated phase boundary curve (magnetic 
field vs temperature). 

The exact solution for the liquid-vapor coexistence curve {pi°y^ vs temperature) of the 
triplet interaction kagome lattice gas (3.1) is given by ^ 

PiT = 1 - i [l - (e""' - 1)^'^'] [1+ < /^o/Ui >L'=o T2<fi>s], < KsjKs < 1, (5.4) 

where < poPi >l*=o, < P >s are the nearest-neighbor pair correlation and spontaneous 
magnetization, respectively, in a standard 5* = | honeycomb Ising model ferromagnet with 
nearest-neighbor pair (dimensionless) interaction parameter K* > and (dimensionless) 
magnetic field L*. It is well known ^ that a necessary and sufficient condition for the 
existence of a phase transition in the 5* = | ferromagnetic {K* > 0) honeycomb Ising model 
is the joint condition L* = and K* > K*, where the critical value K* = | ln(2 + -\/3) = 
0.65847- • •. In (5.4), the exact solutions for < popi >l*=o, and < p >s are known to be 

2 

< /io/ii >L-=o= 3 [coth 2K* + -fKi (k)] , (5.5a) 
with Ki{k,) being the complete elliptic integral of the first kind 

f'7r/2 



and where 



kak)= r {i-n'sm^ey^/He, (5.6a) 

Jo 

= 16^3(1 + z^){l - z)-''{l - z'')-\ (5.6b) 

7 = (1 - z^){z^ - Az + 1) /n\l - z^\{l - zY, (5.6c) 



z = e-"^ . (5.6d) 

The expression (5.4) for pf°^^ can be written completely in the natural /^s-notation of the 
lattice gas by substituting the interaction parameter relation 0] 

K* = - ln(e^^ - 1) at L* = (5.7) 
6 
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into (5.5) and (5.6). Also, the earlier stated critical value K^^ = ln[(2+A/3)^+l] = 3.96992 ■ ■ • 
is determined by substituting the known critical value K* = |ln(2 + a/3) into (5.7). 

The goal of the present section has been attained. Namely, for a given value of the 
single- ion anisotropy parameter a, the longitudinal-transverse coexistence curve {m1°^^ vs 
temperature) of the S = 1 triquadratic Ising model is exactly calculable in m — x space by 
substituting (5.3-7) into (5.2) and replacing whenever appearing by in compliance 
with the mapping relation (3.6b). The resulting coexistence curves display novel behaviors 
as illustrated in Figures 9, 10 and 11. 

In Figure 9, the longitudinal-transverse coexistence curve is plotted for a = —1.5, ex- 
hibiting an asymmetric rounded shape. For increasing temperatures, the transverse (lower) 
branch m™^^ rises more rapidly than the longitudinal (upper) branch mf°'^^ falls, causing 
asymmetry in the rounded shape of the coexistence curve. The resulting curvilinear diame- 
ter of the coexistence region (solid line) increases monotonically, which is more pronounced 
at temperatures closely below the critical temperature (the diameter of the coexistence re- 
gion is defined as the arithmetic mean ^{m^°'^^ + ml°^^) of the longitudinal and transverse 
branches of the coexistence curve or, geometrically, the locus of the midpoints of the vertical 
"tie- lines" spanning the coexistence region). The coexistence curve in Figure 9 is singular 
at its critical point (solid circle). This singular behavior originates within the factor p™^^ of 
the product representation (5.2). Specifically, upon inspecting (5.4), the coexistence curve 
superposes the algebraic branch point (e^/^) and weak energy-type (elne) singularities car- 
ried by < /i >s and < HqUi >l*=o, respectively, resulting in an infinite (vertical) slope at 
its critical point. 

Increasing values of the anisotropy parameter a alter the shape and orientation of the 
coexistence curve. Including the previous case a = —1.5, Figures 10 and 11 exhibit coex- 
istence curves and their curvilinear diameters, respectively, for increasing values of a. To 
illustrate, a = —0.87 reveals an asymmetric rounded coexistence curve (curve 2 in Figure 

10) with a monotonically decreasing curvilinear diameter (curve 2 in Figure 11), and a nar- 
rowing of the coexistence region within the temperature interval 0.5 ^ x < 1, where the 
upper (longitudinal) branch of the coexistence curve falls more rapidly. As argued earlier 
for Figure 9, the coexistence curve 2 in Figure 10 is singular at its critical point (solid circle) 
with the leading singularity again being an algebraic branch point (e^/®). 

As a ascends in value, the m-coordinates of the critical points (solid circles) descend 
in value until vanishing for the value a = a™"^ = —0.81929 ■ ■ ■ (curve 3 in Figures 10 and 

11) . Here the critical point (m > 0) meets its reflective image critical point (m < 0) at 
a confluence point {x,m) = (1,0), thereby changing the topology of the phase diagram. 
Note then that the transverse branches of the coexistence curves entrap a homogeneous 
transverse phase within a tadpole shape region whose positive portion (m > 0) appears in 
Figure 10. Also, the analytical behavior of the coexistence curve near its confluence point 
(x, m) = (1, 0) deserves careful scrutiny since both factors in the product representation (5.2) 
are singular. Sufficiently close to the confluence point, the factor pf°^^ = Pc + o(e^/*) where 
the critical density pc = |(| — v^) = 0.58931 • ■ ■ [0) and using (4.2b), the confluence factor 
[tajohhy^'^ is evaluated along the a = a™"^ phase boundary curve (pbc) as [tanh hy^^ = 
(tanh Q^y^'^ = Q3^^^gi/2_|_Q(^g3/2'j^ (discussed in Section IV. Hence, substituting the above 
expansions into the product representation (5.2), one concludes that, closely approaching its 
confluence point, the a = a^°-^ coexistence curve vanishes as A^e^l"^ (algebraic branch point 
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singularity) where the amplitude Am — PcQ^cA^ — 0.49043 • • •, substituting known values of 

the constants. 

Continuing to increase a beyond its a™"^ value further shortens the length of both the 
coexistence curve (Figure 10) and the corresponding curvilinear diameter (Figure 11) whose 
confluent singularities (open circles) move toward the origin, carrying the same exponent 1/2 
and, suggestively from Figure 10, increasing (decreasing) amplitude for the upper (lower) 
branch of the coexistence curve. Concurrently, the region encasing the homogeneous trans- 
verse phase shrinks in size as it shifts toward the origin. Eventually, as anticipated, phase 
diagrams for the model are nonexistent if a > —2/3. 

The following concepts further elucidate the nature of the confluent phase transitions. 
Since the magnetization order parameter (length of vertical tie-line spanning the coexistence 
region in Figure 10) vanishes continuously at a confluence point, the phase transition is 
a continuous type. One can also argue that the confluence point (open circle) is non- 
critical since it is a non-vertical inflection point of the x = x^onfi isothermal curve in the m 
(magnetization) vs ^ (magnetic fleld) plane. In other words, the initial isothermal magnetic 
susceptibility remains flnite at the confluence point (open circle), contrasting the case of 
a diverging susceptibility which is a thermodynamic hallmark of a magnetic critical point 
(solid circle). The "hybrid" singular point (shaded circle) is also associated with a divergent 
susceptibility as discussed in Appendix B. 

In the phase diagrams of Sections IV and V, we have chosen to accentuate the dipolar 
stimulus ^ and the dipolar thermal response < Si >, both being readily measurable quanti- 
ties. Comparatively, the quadrupolar quantities a and < Sf > are not as easily adjustable or 
accessible experimentally. Theoretically, however, (4.2) directly yields the phase boundary 
curves a(anisotropy) vs a; (temperature) for flxed values of ^(magnetic fleld). Simple geo- 
metrical inspection of Figure 3 reveals that these quadrupolar phase boundary curves are 
devoid of confluent singularities and that each curve ends at a critical point. Also, compos- 
ite differentiation of the mapping relations (3.9) provides the correspondence relationship 
< Sf >= p, leading to a solitary coexistence curve {Sf)'^^ '^ vs a;(temperature), depending 
only upon the temperature in the transition region of Figure 3. 

Lastly, one notes that the single-spin configurational probabilities po, in obvious 
notations, may be represented in terms of the thermal averages < Si >, < Sf > (ordering 
parameters) as p± = |(< Sf > ± < Si >), po = 1— < Sf >. Hence, within the current 
theoretical framework, these conflgurational probabilities can be exactly evaluated along the 
coexistence curves (or across the associated phase boundary curves). Indeed, the resulting 
configurational probabilities confirm earlier descriptions of the S — 1 spin orientations as 
preferentially longitudinal (transverse) in the upper (lower) phases of the diagrams. 

VI. Concluding Remarks 

Exact results in physics are valuable for a variety of reasons. Endeavoring to retain 
only the most essential ingredients of a physical problem, exact solutions of simple model 
systems often provide definite guidance and insights on more realistic and invariably more 
mathematically complex systems. Exact results in tractable models of seemingly different 
physical systems may alert researchers to significant common features of these systems and 
actually emphasize concepts of universality. In addition to their own aesthetic appeal, exact 
results can, of course, serve as standards against which both approximation methods and 
approximate results may be appraised. Also, the underlying mathematical structures of 
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exactly soluble models in statistical physics are rich in content and have led to important 
developments in mathematics. 

In the present theoretical investigations, a S = 1 kagome Ising model was taken to 
have localized triquadratic interactions (J3), single-ion anisotropy {T>), and was placed in 
a uniform magnetic field hz- As anisotropy sufficiently competes with interactions {T> < 
— 2J3/3), the model admits phase transitions and its phase diagrams were determined using 
transformation and mapping methods upon its partition function. The mapping techniques 
significantly simplified the calculations. More particularly, since exact phase diagrams are 
known for condensation in a one-parameter kagome lattice gas model with triplet interactions 
1^, the mapping directly afforded the corresponding exact magnetic phase diagrams of the 
three-parameter S = 1 kagome Ising model. Otherwise, as previously demonstrated with 
the aforementioned lattice gas model, the theory would incorporate a symmetric eight- vertex 
model on the honeycomb lattice in a mediating role and the computations then lengthen 
considerably. 

Special features in the results included the appearance of conffuent singularities causing 
changes in the topology of the phase diagrams, phase boundary curves (magnetic ffeld vs 
temperature) with purely positive or negative slopes as well as intermediate cases show- 
ing nonmonotonicity, and coexistence curves (magnetization vs temperature) with varying 
shapes and orientations, in some instances entrapping a homogeneous phase. The phase dia- 
grams indicate both discontinuous and continuous transitions, the latter being conffuent type 
(open circles), conffuent-critical type (shaded circles), and critical type (solid circles). More 
explicitly, as a function of the temperature, the magnetization order parameter vanishes 
with exponent 1/2 at a conffuent or conffuent-critical singularity, and with exponent 1/8 
at a critical singularity. Also, the magnetic susceptibility diverges at the conffuent-critical 
point with an exponent 7' = 3/4. 

Towards these ends, the phase boundary surface in Figure 2 (temperature as a function 
of the magnetic and anisotropy ffelds) was projected upon the ffeld plane in Figure 3 showing 
the isotherms < T < Tc in the plane to possess crossing points. Such a region of crossing 
points for ffnite ffelds does not occur in the model having biquadratic interactions studied by 
Wu [Q. Hence, some features found in the present phase diagrams may be deemed multi-spin 
interaction effects. As examples, both the purely positive and nonmonotonic slopes of the 
phase boundary curves in Figures 5-7 are connected with the presence of crossing points in 
Figure 3 and their originating triquadratic interactions in the model Hamiltonian (2.1). To 
our knowledge, the exact phase diagrams of the present S = 1 triquadratic Ising model are 
the first obtained for any S = 1 planar Ising model with multi-spin interactions. 
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Appendix A: Crossing points of the a; = isotherm with the x = 1 isotherm 
in the ^ — a plane 

Consider first the ^ > portion of the transition region in Figure 3. As s ^ (T — > 0), 
the leading exponential behaviors in (4.2a) provide the zero temperature isotherm as 
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X 



3a; 



X 



(A.l) 



At criticality, the temperature variable x = 1 {T = Tc) and (4.2a) affords tlie critical 
isotherm as 



1 



pln([(e«3c_i)2/3_ij2cosfi(g3c0) 



Q: 



In 



3c 



4(3 + 2^3) cosli(Q3cO 



(A.2) 



liaving substituted (e^^c _ 1)2/3 _ ^ ^ 2(3 + 2^3) via (4.2b). To locate their crossing point, 
one equates above expressions ag = <^c yielding 

g(€c.<,..+2/3)Q3c ^ 4(3 + 2V3) COsh(Q3cec.o..) 



or 



Across — (2Q3c) ^ In 



2(3 + 2^/3) 



_e2Q3c/3_2(3 + 2y3)_ ' 



(A.3) 



Substituting (4.2b), expression (A.3) gives 

12.92820' 



0.12594- • -In ■ 



1.17806- •• 

which is substituted into (A.l) yielding 

= -2/3 - Across = -0.96837 



0.30171 



(A.4) 



(A.5) 



Similar calculations for the ^ < portion of the transition region in Figure 3 (or by simply 
recognizing the even symmetry of the ^-dependence in (4.2a)) determine the coordinates of 
the crossing points 



±0.30171- ••, tteross = -0.96837- ••. 



(A.6) 



In the vicinity of the crossing points (A.6), countless other crossing points occur among the 
isotherms as shown in Figure 4. 

Appendix B: Divergent susceptibility (7' = 3/4) at the confluent-critical point 

The correspondence relation (5.1) is written 



m = 11 ■ p ■ tanh /i. 



(B.l) 



re-entering the electronic intrinsic magnetic moment jl. Differentiating (B.l) with respect to 
the magnetic field provides an additional correspondence relation between thermodynamic 
response coefficients of the S = 1 triquadratic Ising model and the triplet-interaction kagome 
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lattice gas model. Attention will ultimately be directed toward the "hybrid" singular points 
(shaded circles) in Figure 8 and Figure 10. 

Specifically, the initial isothermal magnetic susceptibility xt of the S = 1 triquadratic 
Ising model is defined (in appropriate units) by 



XT = lim 



dm 



(3ii lim ^ 



da lim 



d{p ■ tanh h) 



3u lim 



dh 

dp didu) , , (9 tanh /i 
■ tanh h + p- 



d{l3p) dh 



dh 



a lim 



dp 
dp 



tanh^ h + Pp/ cosh^ h 



(B.2) 



having used (B.l), the mapping relation (3.6b) and the (dimensionless) field parameter 
h = Pphz- In (B.2), all partial derivatives are performed at constant temperature. 

The grand canonical partition function of the triplet-interaction kagome lattice gas model 
is equivalent (aside from known pre-factors) to the magnetic canonical partition function of 
a standard S = 1/2 honeycomb Ising model with pair interactions and field 0. The thermal 
behaviors of the models exhibit each to have a single critical point with corresponding critical 
exponents being equal. The exponent equivalence will be used shortly regarding the critical 
exponent 7' = 7/4. The isothermal compressibility kt of the triplet-interaction kagome 
lattice gas can be found (in appropriate units) from 




(B.3) 



which in turn is substituted into (B.2) yielding 



XT 



lim 



P^kt tanh^ h + jSp/ cosh^ h 



(B.4) 



the sought relationship between response coefficients xt and kt- 

Consider, in particular, the singular point (shaded circle) in Figure 10 having a confiuence 
of critical (x = 1) points for a = As x — > 1~ (T — > T~) along the a = «^°^ coexistence 

curve, one argues sufficiently near the singular point that kt ^ 00 as e~'^l^ in (B.4), and 
recalls that tanh^ /i — >■ as e along the associated phase boundary curve in Figure 8. Hence 
the relationship (B.4) directly implies that xt ~^ 00 as e~'^l'^ . In effect, the strongly divergent 
(7' = 7/4) susceptibility found at criticality in a standard S* = 1/2 planar Ising model is 
weakened (7' = 3/4) at the confluence of critical points in the S = \ triquadratic Ising model. 
Additionally, the relationship (B.4) verifles that xt remains flnite at the other confluence 
points (open circles) in Figure 10. 
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FIG. 1. The kagome lattice is a two-dimensional periodic array of equilateral triangles and 
regular hexagons. Sites i,j,k are vertices of an elementary triangle. 
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FIG. 2. Phase boundary surface x = x{^,a) where x is the reduced temperature, ^ is the 
reduced magnetic field and a is the reduced anisotropy parameter. Figures 3-8 are projections of 
this surface onto the ^ — a and £, — x planes and will illustrate several unusual features of the surface. 
The reduced parameters ^ = flhz/Js, a = V/J^, x = T/Tc where ksTc/ J3 = Q^^ = 0.25189 • • •. 
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FIG. 3. Isotherms in the ^(magnetic field) — a(anisotropy) plane for several values of 
x (temperature), including the critical isotherm x = 1. The existence of surface undulations in 
FigTirc 2 leads to the crossing of isotherms (boxed region). The reduced parameters ^ = fihz/Ja, 
a = V/J3, X = T/Tc where ksTj J2. = Q^^ = 0.25189 • • •. 
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0.29 0.3 0.31 



FIG. 4. Boxed region in Figure 3 magnified to show details. The reduced parameters 
^ = fihz/Jz, a = P/Js, X = T/Tc where ksTd Jz = Qsc = 0.25189 • • •. 
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0.5 1 

X 

FIG. 5. Projection of Figure 2 onto the ^(magnetic field) — x(temperature) plane for anisotropy 
parameter a = —1.5 and ^ > 0. The phase boundary curve is the locus of discontinuous transitions 
between the longitudinal upper phase and the transverse lower phase, and the solid circle is a critical 
point associated with a continuous transition. The reduced parameters ^ = jlhz/Js, a = "D/J^, 
X = T/Tc where ksTd h = Q^c = 0-25189 • • •. 
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FIG. 6. The ^(magnetic field) — a; (temperature) phase boundary curves for values of the 
anisotropy parameter a where the curves change from being monotonically increasing to monoton- 

ically decreasing, via a non-monotonic region. At across = — -96837 • • •, the a; = isotherm crosses 
the X = 1 isotlierm in the ^ — a plane. The reduced parameters ^ = jlhz/Js, a = V/J^, x = T/Tc 
where ksTc/ h = Qsc = 0.25189 • • •. 
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FIG. 7. The curve a. = across = —0.96837 • • • of Figure 6 magnified to show details of 
the nonmonotonicity. The reduced parameters ^ = jlhz/Js, a = V/J^, x = T/Tc where 
A;Bre/J3 = Q3,' = 0.25189 
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X 

FIG. 8. The ^(magnetic field) — a; (temperature) phase boundary curves for several values of the 
anisotropy parameter a around a^"^ = —0.81929 • • •. The figure extends symmetrically for ^ < 0. 
Curves 1 and 2 end at solid circles representing critical points, while curves 4, 5 and 6 end at open 
circles representing confluent singularities. For curve 3, the critical point meets its reflective image 
critical point at the confluent-critical point (shaded circle). The reduced parameters ^ = jlhz/J^, 
a = P/Jg, X = T/Tc where ksTc/J^ = Q^c = 0.25189 • • •. 
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FIG. 9. The longitudinal-transverse coexistence curve (dotted line) and the corresponding 
curvilinear diameter (solid line) for a = —1.5 in the m(magnetization) — x (temperature) plane. 
The coexistence curve is singular at its critical point (solid circle). The reduced parameters 
m=<: Si>, a = V/Js, X = T/T^ where ksTc/ J2, = Qsc = 0.25189 • • •. 
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FIG. 10. Longitudinal-transverse coexistence curves in the m(magnetization) — a:; (temperature) 
plane for several values of a around a^"^ = —0.81929- • •. Solid circles are critical points, open 
circles are confluence points, and the shaded circle is a confluent-critical point. The flgure extends 
symmetrically for m < 0, so that curves 3, 4 and 5 entrap homogeneous transverse phases within 
tadpole-shaped regions. No such trapped transverse phases exist for curves 1 or 2. The reduced 
parameters m =< >, a = P/J3, x = T/T^ where ksTcl h = Q'^l = 0.25189 • • •. 
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FIG. 11. Curvilinear diameters corresponding to the coexistence curves of Figure 10, showing 
the diameters changing from monotonically increasing to monotonically decreasing behavior and 
ending cither at a critical point (solid circle), or eventually at a confluence (open circle) point or a 
confluent-critical (shaded circle) point. The reduced parameters m =< Si >, a = T)/ J^, x = T/Tc 
where ksTc/ h = Qsc = 0.25189 • • •. The value a™"^ = -0.81929 • • •. 
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